Genome Sequence of Feline Papillomavirus Strain P20 Assembled from Metagenomic Data from the Skin of a House Cat Owner

ABSTRACT A feline papillomavirus genome was assembled from metagenomic sequencing data collected from the skin of a house cat owner. The circular genome of strain P20 is 8,069 bp in length, has a GC content of 54.38%, and displays genome organization typical of feline papillomaviruses. The genome exhibits approximately 75% sequence similarity to other feline papillomavirus genomes.

P apillomaviruses broadly infect mammals, reptiles, and birds (1). Conservation of the major capsid protein L1 suggests that the evolution of papillomaviruses mirrors the phylogeny of their hosts (2). Recent studies have reported on the discovery and diversity of papillomaviruses associated with domestic and wild cats (3). Seven types of feline papillomaviruses, which may or may not produce both skin and oral squamous cell carcinomas in house cats, are currently recognized (4).
We previously studied the diversity and stability of DNA viruses on human skin by sampling three anatomical locations, i.e., left hand, right hand, and scalp, longitudinally over a 6-month period for 43 human subjects (5). Briefly, samples were collected using a tandem dry and wet swab technique. Swabs were then saturated in 1Â phosphate-buffered saline (PBS) and centrifuged at 20,000 Â g. The sample eluent was run through a 0.22-mm filter to remove cellular and bacterial contaminants, and the resulting filtrate was used for viral DNA extraction using the QIAmp ultrasensitive virus kit (Qiagen, Hilden, Germany) and whole-genome amplification (WGA) using the TruePrime WGA kit (Syngen Biotechnology, Inc., Taipei City, Taiwan), each according to the manufacturer's protocol. The resulting DNA was sheared to 600 bp prior to library preparation using the NEBNext Ultra II library preparation kit (New England Biolabs, Ipswich, MA, USA) according to the manufacturer's protocol. The libraries were then sequenced as 150-bp paired-end reads on the HiSeq 2500 platform (Illumina, Inc., San Diego, CA, USA). We identified one participant with a metagenome that consistently showed the presence of a feline papillomavirus, who self-identified as being an owner of a domesticated house cat (5).
To investigate further, 15 metagenome samples from that specific human individual were mapped to a Papillomaviridae reference database compiled from NCBI using BBMap v38.94, with a minimum match of 95% base pair similarity and a k value of 13 (6). Mapped reads were processed using khmer v2.0.0 (7) to remove singletons and overly abundant reads.
We assembled a complete circular feline papillomavirus genome, identified as strain P20, that was 8,069 bp in length and exhibited a GC content of 54.38%. Six open reading frames, including the E6 protein, E7 protein, E1 protein, E2 protein, late protein L2, and major capsid protein L1, which are shared among many animal papillomaviruses (1), were annotated with Prokka v1.14.5 (14) using the embedded viral annotation database and Cenote-Taker 2 v2.1 (15) (Fig. 1) (visualized with SnapGene; GSL Biotech LLC, San Diego, CA). Using BLASTn, our genome assembly showed the greatest similarity (with an average of 70.3% shared amino acids) to feline papillomavirus type 2 (family Papillomaviridae, genus Dyothetapapillomavirus) isolated from the skin of a domestic Maine Coon house cat in 2007 (GenBank accession number NC_038520).
Data availability. All raw sequencing data have been deposited in the NCBI Sequence Read Archive (SRA) under the accession number PRJNA754140. Within that BioProject, the 15 sequencing files used in this study had the following accession numbers: SRS9770471, SRS9770472, SRS9770476, SRS9770495, SRS9770500, SRS9770501, SRS9770587, SRS9770588, SRS9770589, SRS9770614, SRS9770616, SRS9770617, SRS9770642, SRS9770643, and SRS9770644. The complete feline papillomavirus P20 genome has been submitted to NCBI GenBank with the accession number OL310516. All metadata, sequences, annotation files, and scripts used here are publicly available and archived at https://github.com/HerrLab/Graham _2021_feline_papilloma_P20.

ACKNOWLEDGMENTS
This work was completed using the Holland Computing Center of the University of Nebraska, which receives support from the Nebraska Research Initiative. This work was supported by the Department of Justice (grants 2017-IJ-CX-0025, 2019-75-CX-0075, and 2019-R2-CX-0048). This funding agency had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
We thank all of the anonymous participants for their contributions to this study. We also thank members of the Fernando laboratory who assisted with advice on laboratory methods.
We have no competing interests to declare.